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Abstract 

A detailed study of correlated scalars, produced in collisions of nuclei and associated with the 
(7-field fluctuations, (5a) 2 =< a 2 >, at the QCD critical point (critical fluctuations), is performed 
on the basis of a critical event generator (Critical Monte-Carlo) developed in our previous work. 
The aim of this analysis is to reveal suitable observables of critical QCD in the multiparticle 
environment of simulated events and select appropriate signatures of the critical point, associated 
with new and strong effects in nuclear collisions. 
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I. INTRODUCTION 



The existence of a critical point in the phase diagram of QCD, for nonzero baryonic 
density, is of fund— significance for our understanding of strong interactions and so 
its experimental verification has become an issue of high priority For this purpose an 
extensive programme of event-by-event searches for critical fluctuations in the pion sector is 
in progress in experiments with heavy ions from SPS to RHIC energies |2j. In |3| we have 
emphasized, however, that in order to reveal critical density fluctuations in multiparticle 
environment, one has to look for unconventional properties in the momentum distribution of 
reconstructed dipions (7r + 7r~-pairs) [9] with invariant mass just above the two-pion threshold. 
In fact, the QCD critical point, if it exists, communicates with a zero mass scalar field (er- 
field) which at lower temperatures (T < T c ) may reach the two-pion threshold and decay 
in very short time scales owing to the fact that its coupling to the two-pion system is 
strong. Obviously, the fundamental, underlying pattern of a-field fluctuations, built-up 
near the critical point by the universal critical exponents of QCD p|, is phenomenologically 
within reach if and only if the study of correlated sigmas, reconstructed near the two-pion 
threshold, becomes feasible. In the present work we perform a detailed feasibility study of 
the observables related to the detection of the QCD critical point in nuclear collisions. In 
order to proceed we summarize, first, the principles on which the behaviour of a critical 
system of sigmas is based 

QS. 

(a) The geometrical structure of the critical system in transverse space (after integrating 
in rapidity) consists of tr-clusters with a fractal dimension dp = 2 ^7^ leading to a power 
law, < o 2 >~ jxp^+i, for the a-field fluctuations, within each cluster (5 : isotherm critical 
exponent) 

(b) In transverse momentum space the ex-fluctuations obey a power law < a 2 >~ 

2(6-1) 

\p±\ s + 1 leading to observable intermittent behaviour of factorial moments: F 2 (M) ~ 
(M 2 )fer where M 2 is the number of 2D cells j^. 

(c) The density fluctuations of dipions (7r + 7r~-pairs) with invariant mass close to two- 
pion threshold (2m n ) incorporate the sigma- field fluctuations, {5a) 2 ~< a 2 >, at the critical 
point, under the assumption that the sigma mass reaches the two-pion threshold (m CT ~ 2m. K ) 
in a time scale shorter than the relaxation time of critical fluctuations. 

(d) The QCD critical point belongs to the universality class of the 3D Ising system in 
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which 8 « 5. 

On the basis of these principles and the fact that critical clusters in the above universality 
class interact weakly one may construct a Monte-Carlo generator (Critical Monte-Carlo: 
CMC) able to simulate events of critical sigmas, correlated according to the above pre- 
scription 3]. Then, if the corresponding invariant mass distribution is known, we can also 
simulate the decay of the critical sigmas to pions which are experimentally observable. To 
complete the numerical experiment we use the momenta of the charged pions, obtained from 
the decay of the critical sigmas, to form, in an event-by-event basis, neutral dipions (tt + 7t~ 
pairs) and to look for fingerprints of the original critical sigma fluctuations in their momen- 
tum distribution. To check the applicability of our approach we perform also a comparative 
study between the revealed critical correlations-fluctuations in CMC and the corresponding 
behaviour of a conventional Monte-Carlo (HIJING). 

The input parameters of the simulation are the size of the system in rapidity (A) and 
transverse space (R±) as well as the proper time scale (r c ) characteristic for the formation 
of the critical system. In what follows we will use exclusively A = 6 for the rapidity size, 
adapted to the SPS-NA49 conditions (E\, eam ~ 158 GeV/n). The parameters R±, r c can 
then be tuned in order to fit the mean charged pion multiplicity of systems of different size 
(CC, SiSi and PbPb) studied experimentally at this energy p. For each choice of input 
parameters we produce, using the CMC generator, a large set of critical sigma events. In 
order to generate the corresponding pion sector we assign to the sigmas the invariant mass 
distribution p(m cr ) and then we let them decay into pions with a branching ratio 1 : 2 
for neutral to charged. The choice of p(m a ) is determined by the requirement that the 
resulting inclusive neutral dipion invariant mass distribution p(m n + n -) resembles to a large 
extent the corresponding distribution obtained in experiments with heavy ion collisions at 
high energies. In order to perform our analysis in terms of observable quantities we choose 
the mean multiplicity of positive charged pions per event < n n + > as the basic parameter 
(instead of R±, r c ) characterizing different A + A systems. We note that within the CMC 
approach the property < n n + >=< n^- > is exactly fullfiled. We will investigate three 
different cases: low (< n T + >~ 10), intermediate (< >~ 30) and high (< n n + >m 220) 
positive charged pion multiplicity resembling the C + C, Si + Si and Pb + Pb system at 
158 GeV/n respectively. 
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FIG. 1: The critical sigma invariant mass distribution p(m a ) (solid line) as well as the correspond- 
ing neutral dipion distribution / o(m 7r + 7r -) (full circles) for 30000 CMC events with < n n + >= 11.31. 

II. THE CRITICAL SIGMA SECTOR 

First we consider the case of low pion multiplicity. Before going on with the analysis it 
is worth to emphasize that, in general, the determination of the critical sigma sector using 
the observed momenta of charged pions, is a very difficult task. This is due to the absence 
of a characteristic pattern in the inclusive dipion invariant mass distribution attributed to 
the presence of critical sigmas. To illustrate this property we show in Fig. 1 the inclusive 
distribution of the sigma invariant mass calculated before the decay of the sigmas into pions 
(solid line) as well as the corresponding distribution for neutral dipions recontsructed from 
the final pion momenta (full circles) for 30000 CMC events with < n n + >= 11.31 (dataset 
I). Both distributions are equally normalized. Apart from a peak at m^^- m 450 MeV 
due to kinematics there is no other structure in the invariant mass profile. Therefore the 
detection of the critical sigma sector has to go through the, more subtle, study of density 
fluctuations in momentum space. 

In order to reveal the underlying critical fluctuations, at the level of observation, one has 
first to perform factorial moment analysis in small cells of the momentum space [4j. We have 
chosen transverse momenta for this analysis in order to avoid additional assumptions about 
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FIG. 2: The second factorial moment F 2 in transverse momentum space of both the negative pions 
as well as the critical sigmas. We use the same dataset as in Fig. 1. 

;he role of longitudinal rapidity in the description of the statistical mechanics of the system 
jj. Applying factorial moment analysis to the transverse momenta (p x ,Py) of the negative 
pions in the sample of the 30000 critical events we obtain for the second moment a weak 
intermittency effect: F 2 ~ M 2s2 with s 2 ~ 0.077 much smaller than the expected to occur 
in the critical system (s 2 = 2/3). We note with M the number of bins in each momentum 
space component. The corresponding factorial moment for the sigmas, before their decay, 
follows closely the theoretical prediction: s 2 ~ 0.66. This behaviour is displayed in Fig. 2. 

The reason for the suppression of fluctuations in the pionic sector is the kinematical 
distortion of the self-similar pattern formed in the sigma sector due to the sigma-decay. 



The strength of this distortion increases with momentum transfer: Q = J m 2 — Am 2 , and 
becomes negligible near the two-pion threshold [m a « 2m^). Here m 2 = (p n + + p n ~) 2 where 
p w ± are the four momenta of the charged pions produced through the sigma-decay. Thus 
the search for critical fluctuations is based on an accurate reconstruction of the momenta of 
the decaying sigmas. In particular (tt + , 7r~) pairs with invariant mass close to the two-pion 
threshold are the best candidates to carry potentially the geometrical features of the critical 
isoscalars. The great advantage of our approach is that it allows a self-consistency test for 
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the reconstruction of the sigma sector at the level of density fluctuations due to our exact 
knowledge of the sigma momenta and the corresponding fluctuation pattern. 

III. RECONSTRUCTION OF POWER-LAWS 

In practice the reconstruction of the critical momentum fluctuations is performed by 
looking, event by event, for (7r + ,7r _ ) pairs fullfiling the criterion A: 

A = {(n + ,n-)\4ml<(p w+ +p n -) 2 <(2m n + e) 2 } (1) 

with e <C 2m n . The momentum of the corresponding neutral dipion is then obtained as: 
p nn = p 7T + + p^- . In order to ensure that all the available critical sigmas within the above 
kinematical region are recovered in the reconstruction we have to use full pairing forming all 
possible pairs (ti + , 7r~) fullfiling (0) for a given 7r + . However the full pairing introduces as a 
side effect a combinatorial background which has to be treated appropriately. We will come 
back to this point later on. Thus for each value of e we obtain a set of events including dipion 
momenta. We perform factorial moment analysis in transverse momentum space for each 
such dataset. According to the previous discussion, for decreasing values of e the fluctuations 
measured by the intermittency exponent S2 of the corresponding factorial moment should 
increase leading to S2 — > § as e — > + . In Fig. 3a we show the second moment in transverse 
momentum space obtained from data sets of reconstructed dipions for three different values 
of e (5, 50 and 500 MeV). We observe an increase of the slope S2 with decreasing e. However 
this increase has two different origins: (i) the presence of nontrivial fluctuations of dynamical 
origin and (ii) kinematically induced fluctuations through the constraint (|TJ which become 
of the same order as the dynamical ones for e — > 0. It must be noted that the functional 
forms of F2 for the different e values are not exact power-laws and the exponent S2 is an 
effective one. Increasing e the deviations from a power-law description increase too. In 
order to control the kinematical fluctuations we form datasets with mixed events through 
appropriate shuffling of the momenta in the original CMC events. Then we calculate the 
second factorial moments for the datasets consisting of mixed events. We show in Fig. 3b 
the results of this anal ysis . Also in the case of mixed events we observe a clear increase of 
the effective slope [10] for decreasing e. Due to the absence of dynamical fluctuations in 
the mixed events it is natural to conclude that the observed behaviour in F^ 1 ^ is attributed 
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FIG. 3: (a) The second factorial moment F2 in transverse momentum space of reconstructed 
dipions for e = 5 , 50 and 500 MeV (log-log plot). We use the same 30000 CMC events as in 
Fig. 1. In (b) we show for comparison the corresponding factorial moments for datasets consisting 
of mixed events. 

to the kinematical fluctuations. 

These can be removed if we introduce the difference AF q = F q — between the 
factorial moments of real and mixed events (for any order q). For q = 2 the quantity AF2 
represents the correlator of the corresponding dipions. Thus in the difference F2 — F^"^ the 
combinatorial background representing the uncorrelated part of the reconstructed dipions 
is suppressed and the critical fluctuations associated with the correlated part of the dipion 
sector are recovered. We illustrate how this subtraction succeeds in practice in Fig. 4 where 
we present AF2 for e = 5 MeV and < n n + >= 11.31. Although the moments F2 and F^ 
are not exact power-laws, as indicated above, the difference AF2 is significantly better fitted 
by a power-law and the corresponding slope 2 = 0.65(02) is very close to the slope obtained 
from the dataset of the sigmas before their decay (see Fig. 1). 

The fluctuations in the CMC events, as described by AF 2 , are induced by strongly corre- 
lated critical QCD dynamics and are expected to be at least one order of magnitude greater 
than fluctuations originating from conventional hadronic dynamics. To check the validity 
of this statement we have calculated AF2 using 33176 events obtained from the HIJING 
Monte-Carlo generator [7j simulating the SPS C + C system at 158 GeV/n and involving 
only noncritical QCD dynamics. In Figs. 5a-d we show the results of our calculations for 
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FIG. 4: (a) The log-log plot of the correlator AF2 in transverse momentum space of reconstructed 
dipions for e = 5 MeV obtained using 30000 CMC events (the same as in Fig. 1) and the corre- 
sponding mixed events. 

four different values of e. For comparison we plot in the same figure the correlator for the 
30000 CMC events with < n n + >= 11.31. We observe that for all values of e the correlator 
AF 2 for the HIJING system remains flat fluctuating around zero, while for the CMC events 
as e decreases the characteristic, for the critical system, power-law behaviour is setting on 
according to the previous discussions. 



IV. THE CRITICAL INDEX 



It is worth investigating in more detail the role of the parameter e in our analysis. Prac- 
tically the desired limit e — > 0+ is not accessible for several reasons. First of all looking 
at Fig. 1 we observe that in this limit the number of initial critical sigmas as well as the 
number of reconstructed dipions vanishes. We practically need infinite statistics in order 
to extend our analysis in this kinematical region. In addition, in the range of very small e 
values the fluctuations induced by the kinematical constraint (0) become of the same order 
as the dynamical fluctuations due to the critical sigmas. Our entire treatment is based on 
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FIG. 5: The correlator AF2 for the CMC and the HIJING system using four different values of e: 
(a) 5 MeV, (b) 20 MeV, (c) 50 MeV and (d) 120 MeV. 

an accurate cancellation of the kinematical fluctuations in the difference AF2 and this re- 
quires again very high statistics if e — - > 0. In fact, using e values smaller than 4 MeV in the 
analysis of the 30000 CMC events considered so far, the statistical errors increase rapidly 
and prohibit a reliable reveal of the critical fluctuations. On the other hand, if we increase e, 
the combinatorial background due to the reconstructed dipions is increasing rapidly and the 
relative measure of the critical sector goes to zero. Therefore we must search for an optimal 
region of e values to perform our analysis. In order to achieve a more quantitative criterion 
for the determination of the appropriate region of e-values we consider more carefully the 
sector of reconstructed dipions. In the CMC case the reconstructed dipions are divided into 
two subsets: the set of real sigmas and the set of fake sigmas. Let us denote by < n rj(J > e 
the mean number of real sigmas per event with invariant mass in the kinematical window 
for a given value of e. Within our approach we investigate only the part of the sigmas 
which decays into a pair of opposite charged pions. Therefore we have: 

< n r ^+ > e =< n r)7T - > t =< n T>a > t (2) 

where < n r7r +(-) > e is the mean number of positive (negative) charged pions per event pro- 
duced through the decay of the real sigmas with invariant mass in the region (JTJ). Obviously 
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to a good approximation the corresponding number of fake sigmas is given by: 

< n f,a > e ~< n r .^+n r ^- > e — < n r7r + > e (3) 

The first term on the right hand side of eq.Q is dominated by uncorrected (tt + , 7r~) pairs. 
Therefore we can write: < n r ^+n r ^- > e ^< n r>n + > e < n n7r - > t und using the property: 

< n, r 7r + > e =< n r ^- > e we finally obtain: 

< n/ )(T > e ^< n r ^+ >\ - < n r ^+ > e (4) 

A natural constraint ensuring the dominance of the real sigmas over the fake ones is to use 
in our analysis e values fulfilling the condition: 

< nf,a >e < < n r,a >e =^ < n r,n+ >1 ~ < n r,n+ >e < < n r,Tr+ >e (5) 

which simplifies to: < n r 7r + > t < 2. This is fulfilled for example for the value e = 5 MeV 
which we have used so far in our analysis since in this case < n 7T + 7T - > (5 MeV) = 1.12. We 
also observe that < n n + n - > t > < n Tj7T + > e for any value of e, where < n n + n - > t is the 
mean number of dipions per event obtained through the reconstruction in the kinematical 
window (JTJ). The upper bound in < n r 7r + > e corresponds to an upper bound for e. The 
lower bound both for e as well as < n r ^+ > 6 is determined by the statistics according to the 
discussion in the previous paragraph. For practical purposes one applies the stronger bound 

< n n + n - > e < 2 for the estimation of the appropriate region of e- values to be used in the 
data analysis. The great advantage of restricting < ri 7r + 7r - > e instead of e becomes more clear 
when we compare the analysis in CMC datasets simulating A + A processes of different size. 
In this case the dependence of < n^^- > e on e varies from system to system. To illustrate 
this we have produced two additional CMC datasets, with < >= 29.69 (dataset II) 
and < n n + >= 213.96 (dataset III) respectively, each consisting of 30000 events. In Fig. 6 
we show the function < n 1K +^- > t for the CMC datasets I, II and III. The horizontal lines 
correspond to the values < n w + w - > (e) = 5, < n^+ n - > (e) = 1.5 and < n lv + 1T - > (e) = 1.1. 

For increasing size of the system (increasing < 71^+ >) the e values leading to a given 

< n n + n - > e decrease. The complete analysis within our approach for these datasets involves 
the reconstruction of the isoscalar sector as well as the calculation of the corresponding 
correlator AF 2 . However a comparative study between the different datasets in terms of 
factorial moments is possible only for classes of events characterized by almost the same 
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FIG. 6: The function < n^+^- > (e) for the three different CMC datasets described in text. 

multiplicity since in the opposite case artificial fluctuations are induced jjj. Thus, to compare 
the systems of different size, we have to calculate the correlator for fixed multiplicity of 
reconstructed dipions < n^^- > choosing a suitable value of e in each case. As in the case 
of dataset I we use the effective slope (or critical index) <p2 of the correlator (AF2 ~ M 2 ^ 2 ) 
as a measure of the fluctuations in the corresponding dataset. In Fig. 7 we show the results 
of our calculations for the three datasets I, II and III using < n n + n - > e = 5 (full circles), 

< 7V+7T- >e= 1-5 (open stars) and < n^+ n - > e = 1.1 (crosses). We plot 2 as a function 
of the mean number of positive pions per event < n n + > respectively. It is clearly seen in 
Fig. 7 that the critical fluctuations, within the limitations of our analysis due to statistics, 
are at best recovered for CMC systems involving medium to small size nuclei and using 
in the analysis e values fulfilling the constraint < n n + n - > < 2. The critical index 2 
approaches the theoretically expected value §. In fact we see that for < n 7T + n - > t = 1.1 and 

< n n + >= 11.31 we obtain exactly the critical QCD prediction. This is the main result of 
the present work providing us with a useful guide for the search of QCD critical fluctuations 
in relativistic ion collisions. 
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FIG. 7: The slope <p2 as a function of < n w + > for the CMC datasets I, II and III using < n^+^- > 
(e) = 5 (full circles), < n^+^- > (e) = 1.5 (open stars) and < n^+^- > (e) = 1.1 (crosses). 

V. CONCLUSIONS 

In conclusion, we have shown that a set of well prescribed observables (factorial mo- 
ments, correlators, intermittency exponents) associated with the existence of a critical point 
in quark matter, can be established in nuclear collisions. These observables belong to the 
reconstructed isoscalar sector describing massive dipions (it + it~ pairs) near the two-pion 
threshold and their behaviour reveals strong critical effects suggested by u-field fluctuations 
near the critical point. We claim that a search for such a critical behaviour in heavy ion 
experiments is feasible within the framework of a reconstruction procedure of the momenta 
of massive dipions, discussed in this work. The critical effects in this sector although inde- 
pendent of the system size can be at best recovered, through the proposed reconstruction 
algorithm, in collisions of relatively small nuclei. The appropriate kinematical window to 
look for these effects is also determined. Our study has also shown that although it is not 
possible to reveal any conventional sign of the sigma itself in experiments with nuclei at high 
energies, nevertheless its density fluctuations associated with the critical point are observ- 
able and can be measured and studied in a systematic way. Therefore our proposal is to 
study, using the above observables, different processes at the SPS and RHIC with the aim 
to scan a substantial area of the phase diagram, in a systematic search for the QCD critical 
point in collisions of nuclei. 
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